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ABSTRACT 


This  report  describes  a  new  fast  projection  technique  for  adaptive  bearnforming.  It  is 
based  on  a  technique  proposed  by  Yeh  where  Gram-Schmidt  orthogonalization  is  applied  to  N 
rows  of  the  covariance  matrix  to  obtain  the  signal  subspace.  In  Yeh’s  technique,  N  is  equal  to 
L  the  number  of  jammers.  This  technique,  applied  to  bearing  estimation,  requires  a  priori 
knowledge  of  the  number  of  sources.  Even  though  it  gives  go^  results  when  the  numter  of 
sources  is  much  smaller  than  the  number  of  array  elements,  the  performance  of  Yeh's  technique 
degrades  severely  as  the  number  of  jammers  increases.  We  improve  the  technique  by  taking  a 
larger  number  of  tows  in  the  orthogonalization:  we  add  a  thresholding  procedure  to  the  technique 
in  order  to  optimize  N.  We  compared  this  new  technique,  called  Subspace  &.  -ation  without 
Eigenvectors  (SEWE),  with  other  fast  projection  techniques.  SEWE  was  shown  to  give  superior 
performance  in  terms  of  the  Signal-to-Noise-plus-Jammer  Ratio  achievable  with  a  given 
computadonai  load,  except  when  a  very  high  SNJR  is  required. 


r£sum£ 


Ce  report  dferit  une  nouvelle  technique  de  projection  rapide  pour  la  formation 
autoadaptative  des  faisceaux.  Cette  technique  est  basde  sur  une  mdthode  proposde  par  Yeh  ou  une 
orthogonalisation  de  Gram-Schmidt  est  appliqude  h  N  colonnes  de  la  matrice  de  covariance  pour 
obtenir  le  sous-espace  du  signal.  Dans  la  technique  de  Yeh,  le  nombre  de  colonnes,  N,  est  6gai 
au  nombre  de  brouilleurs,  L.  Cette  technique,  utili^  pour  I’estimation  des  directions  des  sources, 
requiert  une  connaissance  a  priori  du  nombre  de  sources.  M6me  si  elle  donne  de  bons  rdsultats 
lorsque  le  nombre  de  brouilleurs  est  tris  faible  par  rapport  au  nombre  d’il^ments,  la  performance 
de  la  technique  de  Yeh  se  degrade  rapidement  lorsque  le  nombre  de  brouilleurs  augmente.  Nous 
ameliorons  la  mdthode  en  lui  ajoutant  un  seuil  ce  qui  permet  d’optimiser  la  valeur  de  N  en 
augmentant  le  nombre  colonnes  utilisdes  pour  faire  Torthogonalisation.  Cette  nouvelle  technique 
est  comparde  avec  d’autres  techniques  de  projection.  On  ddmontre  qu'elle  donne  une  meilleure 
performance  en  terme  du  Rapport  Signal-sur-Bruit-plus-Brouillage  (RSBB)  qu’on  peut  atteindre 
avec  une  charge  de  calcul  doimte,  sauf  dans  le  cas  de  RSBB  trhs  61ev6. 
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EXECTJnVE  SUMMARY 


Jamming  of  radar  systems  via  the  antenna  sidelobes  or  the  antenna  mainbeom  con 
seriously  degrade  radar  performance.  Thus  jamming  can  provide  a  crucial  advantage  to  enemy 
forces  if  radar  electronic  counter-counter-measures  (ECCM)  are  not  effective. 

Adaptive  antenna  nulling  provides  an  effective  ECCM  but  at  the  cost  of  increased 
equipment  complexity  and  capability.  A  major  cost  factor  is  the  need  for  a  very  high  speed,  real- 
time,  computadonai  capability.  The  computational  capability  required  is  so  high  it  limits  the 
application  of  this  effective  ECCM. 

This  report  decribes  and  evaluates  a  computationally  efficient  algorithm  that  reduces  the 
requirements  for  high  speed  computation  by  10  times  in  some  cases. 

The  new  method,  called  Subspace  Estimation  Without  Eigenvectors  (SEWE),  uses  the 
concept  of  signal  and  noise  subspace.  Here  the  noise  subspace  is  the  jamming  space.  The 
algorithm  calculates  a  projection  matrix  which  represents  a  subspace  orthogonal  to  the  jammer 
subspace.  The  projection  matrix  is  applied  to  a  weight  vector  that  steers  the  array  in  a  particular 
direction.  When  the  resulting  adapted  weight  vector  is  applied  to  superimposed  signal  and 
jamming,  the  jamming  signals  are  greatly  reduced  while  the  loss  of  signal  is  small.  This  occurs 
because  the  a^pted  weight  vector  is  approximately  orthogonal  to  the  jamming  signals. 


We  compared  the  SEWE  technique  with  other  fast  projection  techniques.  SEWE  was 
shown  to  give  superior  performance  in  terms  of  the  achievable  Signal-to-Noise-plus-Jammer 
Ratio  (SNJR)  for  a  given  computational  load,  except  for  the  case  where  very  high  SNJR  is 
required. 
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SUBSPACE  ESTIMATION  WITHOUT  EIGENVECTORS  FOR  ADAPTIVE 

BEAMFORMING 


1.0  Introduction 

*  Array  signal  processing  techniques  based  on  eigenvector  projection  ore  highly  effective 
for  both  interference  suppression  and  for  estimation  of  signal  angle-of-anrival.  Projection  methods 
employ  the  concept  of  signal  and  noise  subspaces.  In  the  case  of  jamming  signals,  the  noise 

*  space  represents  the  jamming  space.  Jamming  is  cancelled  by  projecting  the  array  dau  vectors 
into  a  space  orthogonal  to  the  noise  space,  ie  the  signal  space.  In  the  case  of  angle-of-arhval 
estimation,  the  array  data  vectors  ate  projected  into  a  space  orthogonal  to  the  signal  space.  These 
techniques  estimate  the  noise  and  signal  subspace  with  an  error  that  decreases  when  an  increasing 
number  of  sample  vectors  is  used  in  the  calculation  of  the  covariance  matrix,  /^.veraging  is 
exploited  by  incorporating  many  data  vectors  into  an  estimated  covariance  matrix  which 
approaches  the  true  covariance  matrix  asymptotically  as  the  observation  period  increases. 

BUhiing  n>21  proposed  an  eigenvector-based  projection  technique  for  adaptive 
beamforming;  the  eigenvectors  corresponding  to  the  strongest  eigenvalues  of  the  covariance 
matrix  are  used  to  calculate  a  projecdon  matrix  orthogonal  to  the  jammer  subspace.  This 
projection  matrix  is  applied  to  the  array  sreering  vector  which  is  projected  into  a  subspace 
orthogonal  to  the  interference.  The  resulting  adapted  steering  weight  vector  is,  ideally, 
orthogonal  to  the  jamming  signals.  Variations  of  the  eigenvector-based  adaptive  nulling 
algorithms  have  been  studied  by  Gabriel  [3].  He  sho^ved  that  the  sidelobes  of  low  sidelobe 
antennas  were  less  perturbed  when  eigenvector-based  projection  techniques  were  used  as  an 
alternative  to  the  Sample-Matrix-Inversion  (SMI)  technique  of  Reed  et  al  [4].  Because  the  signal 

«  or  interference  subspace  is  estimated  with  a  fmite  sample  size,  the  cancellation  of  the  jamming 

is  imperfect. 

*  Projection  techniques  based  on  eigenanalysis  are  computationally  demanding  particularly 
for  arrays  with  a  large  number  of  elements.  This  limits  the  applicability  of  these  techniques  where 
high-speed  real-time  processing  is  required.  Fast  projection  techniques  based  on  Data  Vector 
Orthogonalization  (DVO)  [S,6,7,8]  avoid  calculation-intensive  eigenanalysis,  but  at  a  cost  in 
performance.  The  loss  in  perfoimance  comes  from  the  inability  of  these  techniques  to  incorporate 
averaging  to  better  estimate  the  signal  or  interference  subspace.  In  contrast,  eigen-based  methods 
become  progressively  more  accurate  as  the  number  of  data  vectors  incorporated  in  the  covariance 
matrix  increases.  In  the  most  advanced  version  of  DVO,  Nickel  [8]  introduces  a  double-threshold 
procedure  that  improves  performance;  an  internal  threshold  rejects  data  vectors  that  do  not  add 
sufficient  information  about  the  signal  or  interference  subspace,  while  an  external  threshold  keeps 
the  procedure  going  until  a  perfoimance  criterion  is  satisfied.  Toulgoat  and  Turner  [9,10,11] 
*'ave  developed  a  technique  called  Data  Vector  Selection  and  Orthogonalization  (DVSO)  which 
provides  a  better  estimate  of  the  interference  subspace  than  does  DVO,  albeit  at  the  price  of  a 
higher  computational  load.  The  required  load  is,  however,  still  much  less  than  that  of  the 
eigenvector  techniques.  In  addition,  Toulgoat  and  Turner  [9,10,11]  have  used  DVSO  with 

.  averaging  procedures  to  develop  two  algorithms  that  provide  a  better  trade  off  in  performance 
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versus  comp  itauonal  load.  The  first  technique,  caJJed  Data  Vector  Selection.  OrthogonaJizatioa 
and  Weight  vVERaging  (DVSO-WAVER).  averages  repeated  computations  of  the  adapted  weight 
vectors  to  '^ive  a  final  adapted  weight  vector  for  nulling  the  jamming.  The  second,  called  DVSO- 
COVAR,  uses  a  data  selection  process  for  the  choice  of  data  vectors  used  in  the  formation  of  the 
covariance  matrix.  Gram-Schmidt  (GS)  orthogonalizalion  is  then  applied  to  the  columns  of  the 
covariance  matrix  to  obtain  a  good  estimate  of  the  signal  subspace.  These  two  methods  provide 
essentially  the  same  Signal-to-Noise-plus-Jammer  Ratio  (SNJR)  as  the  SMI  method  but  require 
far  fewer  computations. 

The  data*vector  selection  technique  used  in  DVSO-COVAR  is  similar  to  a  technique  that 
Reilly  and  Law  [12]  used  for  direction  finding.  However,  [12]  applies  the  selection  process  to 
the  permutation  and  selection  of  columns  of  a  covariance  matrix.  In  contrast,  in  DVSO-COVAR, 
the  data  selection  is  applied  prior  to  the  calculation  of  the  covariance  matrix,  in  order  to  choose 
the  data  vectors  to  be  used  in  the  formation  of  the  covariance  matrix.  A  procedure  equivalent  to 
that  of  [12]  is  then  applied  for  selecting  a  subset  of  the  columns  vectors  of  the  covariance  matrix. 

Yeh  [13,14]  proposed  another  technique  based  on  Gram-Schmidt  orthogonalizalion.  This 
technique  was  first  applied  to  bearing  estimation;  the  projection  matrix  is  computed  by  using  any 
L  rows  of  the  covariance  matrix  where  L  is  the  number  of  jammers.  It  is  computationally  more 
efficient  than  projection  techniques  using  eigenvectors  especially  for  the  case  where  the  number 
of  jammers  is  much  smaller  than  the  number  of  array  elements.  In  diis  report,  the  Yeh  technique 
is  applied  to  adaptive  beamforming  and  shown  to  give  good  results  when  the  number  of  jammers 
is  much  smaller  than  the  number  of  array  elements.  We  propose  an  improved  version  of  the  Yeh 
technique  incorporating  a  threshold  procedure,  called  Subspace  Estimation  Without  Eigenvectors 
(SEWE).  This  new  algorithm  is  compared  with  eigenvector  decomposition.  DVSO  and  DVSO- 
COVAR  techniques. 


2.0  Signal  and  jammer  characteristics 

Radar  systems  usually  transmit  pulsed  signals  that  are  narrowband  with  respect  to  the 
carrier  frequency.  Since  the  targets  of  interest  are  in  the  far  field,  the  received  signals  are  plane 
waves.  When  no  reflection  is  received  from  a  target,  the  data  vector,  x(tn),  is  defined  by: 

L 

X(ta)  -  •  A  J„  *  IXj, 

where  L  is  the  number  of  jammers,  jB=[ji(U.....jic(U]  jammer  vector  with  j^t,)  being 
complex  gaussian  amplitude  of  the  i*  jammer  at  time  t„,  ii„  is  the  receiver  noise  vector  with 
mutually  independent  components.  Its  power  is  E{n,*}=K(J^  where  K  is  the  number  of  array 
elements.  The  sampling  rate  is  selected  such  that  the  quantities  ji(t,)  and  ji(t,^,|)  are  statistically 

independent.  The  matrix  A  represents  the  directions  of  the  jammers  with  A=[a(0,) . aOic)].  The 

quantity  aOj)  is  a  deterministic  vector  representing  the  direction  of  arrival  of  the  i*  jammer. 
Without  any  loss  of  generality,  we  assume  a  uniform  linear  array  with  equally  spaced  elements. 
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In  that  case,  a(d,)  is  defined  as 

*(6^)  *■-[!,  exp (ji^ sin (0^)) _ _ exp  ( j  sin (6,)  )  ]  (D 

where  d  is  the  inter-element  spacing.  6i  is  the  direction  of  arrival  of  the  i'‘‘  jammer,  X  is  the 
jammer  wavelength. 

When  projection  techniques  are  applied  to  radar  systems  for  adaptive  nulling  of  jammers, 
it  is  important  t^  sample  vectors  be  taken  in  the  absence  of  signal  and  clutter.  This  is  a 
standard  radar  problem;  one  solution  is  to  sample  at  time  intervals  corresponding  to  ranges  where 
it  is  Iciiown  that  there  are  no  targets  or  clutter  returns. 


3.0  Performance  criterion 

The  ultimate  measure  of  performance  of  an  adaptive  array  is  the  Signal-to-Noise-plus- 
Januner  Ratio  (SNJR).  We  attempt  to  maximize  this  quantity  by  using  an  adapted  weight  vector, 
w.,  that  discriminates  against  jamming  while  causing  a  minimal  loss  of  useful  signal.  The  adapted 
weight  vector,  w,,  is  calculated  from  a  steering  vector.  w„  which  specifies  the  direction  for  which 
the  SNJR  is  to  be  maximized: 


•r.  -  a  r  ir, 

where  T  is  a  transformation  matrix  and  a  an  arbitrary  constant.  The  optimum  solution  for 
maximisation  of  the  SNJR  is  given  by  T=4>''  where  O  is  the  true  covariance  matrix.  For  an 
infinite  Jammer-tc- Noise  Ratio  (JNR)  in  each  channel  the  transformation  can  be  written  as 

T  -  I-A{A"A)-^A‘’ 


where  A  is  the  matrix  of  the  directions  of  the  jarruners  given  by  (1).  The  transformation  matrix 
T,  projects  the  steering  vector  into  a  subspace  orthogonal  to  the  subspace  which  is  spanned  by 
the  directional  vectors  of  the  jammers  (jammer  subspace).  Instead  of  using  the  duection  vectors 
of  the  jarruners  to  build  up  the  projection  matrix,  T,  one  can  use  any  L  basis  vectors  which  span 
this  subspace.  The  projection  matrix,  T,  is  written  as: 

r  -  r  -  W*  (2) 

where  Vs[v„...,vJ  is  the  matrix  of  the  L  basis  vectors.  When  sensor  noise  is  present  (especially 
for  low  JI^),  we  may  need  more  than  L  basis  vectors  to  better  define  the  jarruner  subspace. 

With  this  adapted  weight  vector  the  SNJR  is  given  by: 
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SNJR{9) 


g  { lt>i[  a(e)P  )  _  la(9)*ra(e)P 

g  {  wj  g  •(0)"  r<^  r  «(0)  k 

where  a(0)  is  the  steering  vector  with  0  being  the  steering  direction.  The  quantity  tr/K  is  the 
normalization  factor. 

We  calculate  a  two-dimensional  SNJR  to  evaluate  the  performance  in  all  directions.  It 
is  defined  as  follows: 


2 

SNJR^  -  Y,  SNJlt(B) 


An  area  of  one  beamwidth  around  each  jammer  is  excluded  since  the  SNJR  degrades  severely 
at  those  locations.  The  beamwidth,  BW.  is  defined  as  follows: 


where  L  is  the  length  of  qieiture;  L  a  K  X/l, 


3.1  Test  case 

The  array  has  64  elements  with  half-wavelength  spacing  between  each  element. 
Performance  is  evaluated  for  three  different  scenarios:  scenario  #1  with  La3  jammers  located  at 
0,8-33.4”,  02S-12.7”  and  038-8.6”;  scenario  #2  with  LsS  jammers  located  at  0,8.41.3”,  02=- 
24.2”,  838-9.2”,  04=-6.3”  and  9j=-3.4”;  scenario  #3  with  LslO  jammers  located  at  0,=-45.2”,  02=- 
41.3”,  038-37.6”,  04=-24.2”,  0j=-18.l”.  0s=-l5.1”,  07=-12.1”,  0,=-9.2”,  0^-6.3”  and  0,o=-3.4”. 
We  use  two  different  values  for  the  JNR:  26.9  dB  and  40  dB.  The  noise  power,  o^,  is  set  at  1. 


4.0  Yeh-Projcction  technique 
4.1  Description  of  the  technique 

The  Yeh-projection  technique  [13.14]  calculates  the  projection  matrix  using  any  L 
columns  of  the  covariance  matrix,  where  L  is  the  number  of  jammers.  The  sample  covariance 
matrix  is  defined  as  follows: 


R 


. . .  ,r,f] 


(3) 
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where  M  is  the  number  of  samples  used  in  (he  formation  of  the  cov  lu.  .cc  matrix  and  r,  is  the 
J’'*  column  of  the  covariance  matrix.  The  jammer  subspace.  /,  is  calculated  by  an 
erthonormalization  of  the  L  columns  of  the  covanance  matrix,  R.  that  have  been  cstimaiod.  The 
projection  matrix  is  then  calculated  with  equation  (2). 

Yeh  { 1 3, 1 4]  used  this  projection  matrix  to  evaluate  a  spectral  estimator  and  then  determine 
the  positions  of  the  sources.  It  gives  poor  results  when  the  Signal-to-Noise  Ratio  is  low  [13]. 
Moreover,  the  direction  finding  application  requires  the  technique  to  be  able  to  discriminate 
between  two  sources  which  are  very  close  together.  In  contrast,  adaptive  beamforming  has  only 
to  cancel  jamming  signals  which  are  very  strong  with  respect  to  the  receiver  noise  power, 
furthermore,  there  is  no  need  for  two  different  nulls  if  the  signals  are  very  close  together.  For 
those  reasons,  the  projection  technique  of  Yeh  is  well  suited  for  jammer  nulling.  The  projection 
mauix  is  used  to  calculate  an  adapted  weight  vector,  w,,  from  a  steering  weight  vector,  w,. 


w.  =  Tw, 


where  the  steering  weight  vector.  w„  is  given  by 

tf,  -  [1,  exp(2-ji  d  sin(6,) ),...,  exp(2-^  {K-D  i  d  sin(0,)  ] 


Besides  its  good  performance  in  terms  of  SNJR,  Yen’s  technique  has  a  small 
computational  load  compared  with  other  adaptive  beamforming  algorithms  such  as  Biihnng 
technique,  SMI  or  the  DVSO-COVAR  technique.  The  computational  load  for  Ych's  technique 
is  given  in  terms  of  the  number  of  complex  multiplications  (c.m.)  required  to  find  the  adapted 
weight  vector 

MKL*L^K 

where  M  is  the  number  of  samples  used  in  the  formation  of  the  covariance  matrix.  K  is  the 
number  of  array  elements  and  L  is  the  number  of  jammzrs.  Note  that  only  L  columns  of  the 
covariance  matrix  have  to  be  calculated,  which  reduces  the  number  of  complex  multiplications 
to  MKL  (c.m.)  compared  with  M(K-»*l)K/2  (c.m.)  for  the  calculation  of  the  lower  triangle  of  the 
covariance  matrix.  If  L«K,  this  greatly  reduces  the  computational  load  incurred  in  the 
calculation  of  the  covariance  matrix. 

The  next  section  presents  simulation  results  for  the  Yeh-projection  technique  as  applied 
to  adaptive  beamforming. 
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4^  Simulation  resuit5*Yeh’s  technique 

In  this  section.  Ych’s  technique  is  compared  with  conventional  Founcr  bearnfoimin^  with 
Taylor  antenna  tapering  to  control  the  sid.:lobes.  The  comparison  is  earned  out  by  scimninji;  the 
array  over  the  enure  range  of  angles.  0.  i.c.  for  -15  sini8)  51.  As  the  scan  angle  passes  the 
jamming  signal  there  is  a  very  large  decrease  in  SNJR;  however  the  region  of  decreased  SNJR 
is.  in  general,  very  much  narrower  fo<  'he  Yeh  technique. 

Figures  1,  2.  3  show  the  no"  '.£ed  SNJR  versus  sin(0)  for  scenarios  1  tLsS),  2  (L=5) 
and  3  (LslO)  for  both  the  Yeh  techu...,t:e  and  Fourier  beamforming  w;th  Taylor  weights.  For  all 
of  the  scenarios,  the  JNR  is  26.9  df  For  the  small  ratio  L/K  in  scenario  1.  the  SNJR  curve  for 
the  Yeh  technique  approaches  the  o:  i'um  value  of  0  dB  except  .'or  areas  around  the  januner. 
The  SNJR  for  Fourier  beamforming  wi*h  Taylor  weighting  is  severely  degraded  everywhere.  As 
the  number  of  Jammers  increases  or  as  they  get  closer  together,  the  performance  of  the  Yeh- 
projection  degrades  progressively  as  seen  in  Ingure  3.  Here  the  number  of  jammers  has  increased 
to  ten  with  a  cluster  of  six  closely  spaced.  In  Figure  3,  the  SNJR  for  the  Yeh  technique  is  8  dB 
below  SNJR  maximum. 

The  performance  of  the  Yeh  technique  can  be  improved  by  using  a  larger  number  of 
columns  in  the  orthonormalization.  To  evaluate  the  number  of  columns  required,  the  normalized 
SNJRjo  is  plotted  versus  the  number  of  columns,  N,  used  in  the  orthonormalization.  Figure  4 
gives  SNJR.0  versus  N  for  scenarios  1  (L»3),  2  (L=5)  and  3  (L=10)  respectively  with  JNR=26.9 
dB.  In  the  first  scenario  where  L=3,  the  optimum  value  for  N,  No^,  equals  the  number  of 
jammers.  L.  For  L=5.  the  optimum  value  for  N  stands  at  Nopi=8.  while  the  value  N^=19  gives 
the  maximum  SNJR-d  when  L=10.  We  observe,  in  general,  that  SNJRy)  is  very  poor  for  N<L. 
and  that  the  SNJR^)  decreases  slowly  as  a  function  of  N  for  N>N^.  Figure  5  gives  the  same 
results  for  JNR=40  dB.  In  this  case,  tends  to  be  slightly  lower  as  expected. 

Figures  6,  7  and  8  display  SNJR  versus  sin(0)  curves  for  N=L  and  NsN^^  for  the  same 
scenarios  as  above  with  the  JNR=26.9  dB.  As  expected,  using  N=N^  instead  of  N=L  greatly 
improves  the  performance  of  the  Yeh  technique  for  scenarios  2  and  3  while  leaving  the 
performance  unchanged  for  scenario  1.  Yeh’s  technique  can  then  be  modified  by  using 
columns  instead  of  L  for  the  Gram-Schmidt  orthogonalization.  Since  the  quantity  N^^  increases 
with  the  number  of  jammers,  it  is  not  possible  to  find  a  single  value  for  N„^  convenient  for  every 
scenario.  In  order  to  determine  N^  for  each  scenario,  we  have  devised  a  thresholding  procedure 
for  Yeh’s  technique.  This  procedure  allows  the  number  of  sample  vectors,  N,  to  increase  beyond 
L  to  approximately,  N^,  eliminating  the  requirement  for  a  priori  knowledge  of  the  number  of 
jammers  and  their  disposition  and  power.  In  the  next  sccL  '^  present  a  new  version  of  Ych’s 
technique,  called  Subspacc  Estimauon  Without  Eigenvectors  (SEWE). 
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sin  (9) 

Figure  1  Normalized  SNJR  versus  sin(8)  for  L=3.  K=64  and  JNR=26.9  dB. 


sln(0) 


Nonnalized  SNJR  versus  sin(6)  for  L=S,  K=64  and  JNR=26.9  dB. 
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Number  of  basis  vectors,  N 

Figure  5  Nonnalued  SNJRjo  versus  N  for  Yeh  with  K=64  and  JNR340  dB. 


sin(0) 

Figure  6  Normalized  SNJR  versus  sin(0)  for  L=3,  K=64  and  JNR=26.9  dB. 
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Normalized  SNJR(dB) 


5.0  Subspace  Esdnution  Without  Eigenvectors  (SEWE) 

The  SEWE  technique  trades  off  a  small  increase  in  the  compuutionai  load  for  an 
improved  SNJR.  It  is  shown  that  the  addition  of  a  thresholding  procedure  to  the 
orthonormaiization  procedure  considerably  improves  the  performance  of  Ych’s  technique  at  the 
cost  of  a  moderate  increase  in  computational  load.  The  thresholding  procedure  has  two 
advantages:  firstly,  a  priori  knowledge  of  the  number  of  jammers  is  not  required;  secondly,  the 
procedure  leads  to  selection  of  values  for  N  close  to  Ng,|. 

We  investigate  two  different  statistics.  The  first  statistic  is  defined  as  the  power  of  the 
residual  basis  vector  after  orthogonalization,  T,  [3,6]: 

^1.*  -  luj*  (^> 

where  u,  is  the  residual  basis  vector  after  orthogonalization  at  the  m"*  step.  This  is  compared  to 
a  threshold.  A,,  which  is  proportional  to  the  noise  power.  The  procedure  terminates  when  the 
power  of  the  residual  basis  vector  is  approximately  equal  to  the  noise  power.  The  second  statistic 
T]j,  is  based  on  the  volume  of  the  generalimd  pai^elepiped  defined  by  the  vertices  of  the 
vectors,  [r|,...jrJ  where  Ti  is  the  i**  column  vectm  of  the  covariance  matrix  [8].  Noting  that  u, 
is  the  i^  orthogonal  basis  vector  generated  from  the  column  vectors,  r„  up  to  Tj,  the  volume 
of  the  parallelepiped  is  given  by: 


It  is  also  noted  that  tiP  ^  kij'.  A  test  statistic  based  on  the  square  of  the  volume  is  given  by 


_  ■  lu/ 

i-1  b/ 


(5) 


We  observe  that  Tjj,  S  1,  with  when  the  r,’s  and  u,’s  are  colinear.  Secondly,  because  the 
u,’s  are  generated  from  the  r,’s  by  the  Gram'Schmidt  process,  Tjj.  is  independent  of  any  overall 
scale  far  tors  or  gain  factors  i.e.  multiplying  all  of  the  data  vectors  by  a  constant  wiU  not  change 
T]^.  In  general  as  m  increases,  a  point  will  be  reached  where  the  jamming  space  is  well 
rejnesented  by  the  Nsm  vectors,  u,  to  0^.  T2jb  decreases  rapidly  as  m  increases  beyond  N.  The 
optimum  value  for  N,  N^p^  is  selected  when  the  quantity,  T2  which  is  the  length  of  the  edge 
of  a  cube  of  volume  equal  to  CT])*^  falls  below  a  threshold,  A2.  Since  the  statistic  Tj  is  always 
sixudler  than  1,  A,  is  chosen  accordingly. 

The  details  of  the  new  algorithm  follow.  As  in  Yeh’s  technique,  the  covariance  matrix  is 
calculated  one  column  at  a  time: 
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(6) 


M 

where  is  the  element  of  the  m'*'  column  of  the  covariance  matrix. 

# 

The  algorithm  is  presented  here  with  a  statistic  T.  and  a  threshold  A,  which  are  equal  to 
either  T|^  and  A,,  or  Tij,  and  A}  depending  on  the  statistics  investigated. 

The  technique  is  in^ilemented  as  follows: 

1)  Set  n^l  and  calculate  r^  for  (si  to  K. 

2)  u,  s  r.  »  [r,.  rj.  ...  r^J^ 

3)  Calculate  a  statistic  T.  according  to  equation  (4)  or  (S) 

4) 

IfT.>A 

else 

T,i  3  0;  the  procedure  terminates  indicating  no  jammers  are  present 

5)  Set  mstiH-l 

6)  Calculate  r^  for  (si  to  K  according  to  (6).  t 

7) 

8)  Calculate  the  statistics  according  to  (4)  or  (5). 

9) 

If7,>A 

else 

=  0;  the  procedure  terminates  with  N  basis  vectors  where  N=m-1. 

Calculate  w,  as. 
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•f,  -  ir,  -  ^  (v^iir,) 


i-l 


The  choice  of  a  value  for  A  depends  on  the  choice  of  the  statistics.  If  T„  s  T,  A  is  chosen 
propoitional  to  noise  power.  If  «  T2jb«  Aj  is  related  to  a  length.  The  addition  of  a  threshold 
in  the  procedure  increases  the  computational  load.  The  new  computational  load  is  given  by: 

€•  Nl  (M  ~  .75)1?  ♦  ^  *  (lf^.5)] 

M  m 


where  N  depends  on  the  threshold  value.  As  the  threshold  A  decreases,  the  computational  load 
increases. 


5.1  Simulation  results  •  SEWE  technique 

We  first  determine  through  simulation  the  optimum  values  for  A,  and  A2.  We  also  analyze 
the  effect  of  the  number  of  jammers  and  of  the  Jammer>to-Noise  Ratio,  JNR,  on  the  selection  of 
the  threshold.  Figure  9  gives  normalized  SNJRjo  versus  A|  for  scenarios  1,  2  and  3  with 
JNR=26.9  dB.  AH  three  scenarios  have  a  maximum  SNJRjo  for  a  value  of  A,^  between  .65  mW 
and  3.2  mW.  As  A,  gets  larger  than  A,^.  the  SNJRjo  decreases  very  quickly  for  scenarios  2  and 
3.  This  region  corresponds  to  a  small  number  of  basis  vectors  (Same  area  as  N<L  in  section  4.2). 
On  the  other  hand,  as  A,  gets  smaller  than  A,^  the  SNJRa>  gently  rolls  off  (Case  where  N>No^ 
in  section  4.2).  Setting  the  threshold’s  value  to  1.3  mW  gives  performance  near  optimum  SNJRjj, 
(±0.1  dB). 

SNJRjo  versus  A|  is  shown  in  Figure  10  for  the  same  scenarios  as  above  with  JNR=40 
dB.  In  this  case,  a  value  of  1.3  mW  for  A,  gives  the  near  maximum  value  (±0.15  dB)  for  the 
three  scenarios.  The  value  of  the  threshold  is  only  slightly  affected  by  the  number  of  Jammers 
and  the  JNR.  However,  the  threshold  depends  on  the  number  of  array  elements,  K,  and  is  also 
propoitional  to  the  squared  power  of  the  noise  power,  p,. 

We  then  investigate  the  second  statistic  T2  and  the  threshold.  Ax-  Figures  1 1  and  12  give 
SNJRjb  versus  A,  for  scenarios  1.  2  and  3  with  JNR=26.9  md  40  dB  respectively.  We  observe 
basically  the  same  behavior  as  for  A|.  A  value  of  0.005  for  A^  gives  near  maximum  SNJRjo  (± 
0. 15  dB)  for  SEWE  technique  for  the  three  scenarios  with  the  two  different  JNR’s.  This  threshold 
is  independent  of  the  noise  power,  a  very  practical  advantage. 

Figures  13,  14  and  15  give  SNJR2D  versus  sin(0)  for  scenarios  1,  2  and  3  respectively 
when  using  Tj.  Results  for  conventional  beamforming  and  the  Yeh-projection  technique  are  also 
shown  on  the  same  graph.  It  is  clear  from  these  results  that  adding  a  thresholding  procedure  to 
Yeh’s  technique  improves  the  performance.  The  next  section  compares  the  SEWE  technique  with 
DVSO-COVAR  and  an  eigenvector-based  projection  technique. 
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Normalized  SNJR^(d8)  I  Normalized  SNJR^(dB) 


(mW) 

Normalized  SNJRjo  versus  A|  for  SEWE  with  K=64  and  JNR=26.9  dB 


(mW) 


Figure  10  Normalized  SNJRjo  versus  A,  for  SEWE  with  Ks64  and  JNRs40  dB. 
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Normafized  SNJR  (dB) 


Figure  11  Nonnalized  SNJRjq  versus  d,  for  SEWE  with  Ks64  and  JNR326.9  dB. 


A2 

Figure  12  Nonnalized  SNJR]],  versus  Aj  for  SEWE  with  Ks64  and  JNRs40  dB. 
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sin(0) 

Nonnalized  SNJR  versus  sin(6)  for  SEWE  with  K=64  and  JNRs40  dB. 


sin(0) 

Figure  14  Nonnalized  SNJR  versus  sin(6)  for  SEWE  with  L=S,  K=64  and  JNR=26.9  dB 
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Figure  IS  Normalized  SNjR  versus  sm(6)  for  SEWE  with  LalO,  K=64  and  JNR=26.9  dB 


6.0  Comparison  with  other  projection  techniques 

The  SEWE  technique  is  compared  with  three  projection  techniques;  DVSO,  DVSO- 
COVAR  and  an  eigen-decomposition  based  projection  method.  A  brief  description  is  given  of 
1  the  projection  techniques.  This  is  followed  by  simulation  results  illustrating  the  comparative 

performance  of  the  various  techniques. 


6.1.1  DVSO  technique 

The  basis  of  the  Data  Vector  Selection  and  Orthogonalization  (DVSO)  [9,  10,  11] 
technique  is  the  selection  of  the  best  M  vectors  from  a  larger  group  of  N  vectors  which  are 
sampled  when  the  transmitting  signal  is  absent  and  stored  before  the  process  begins.  The  M 
vectors  are  selected  on  the  basis  that  they  yield  the  highest  SNJR  as  compared  to  any  other  set 
of  M  vectors  selected  from  the  set  of  N.  The  selection  process,  the  Gram-Semidt 
orthogonalization  and  the  computation  of  the  adapted  weight,  w„  are  combined  in  a  single 
procedure.  The  computational  load  is  given  by  [16]: 

Cbvso  -  J:(0 . 5+0 . 5N*2M*MN*M^)  +0  .  SMN 
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6.1J!  DVSO-COVAR  technique 


DVSO-COVAR  [11]  employs  the  DVSO  method  in  calculating  the  covariance  matrix  and 
processing  its  columns  to  obtain  an  adapted  weight  vector.  The  major  feature  of  DVSO-COVAR 
is  its  selection  process:  the  best  data  vectors  are  chosen  to  form  the  covariance  matrix.  We  are 
thus  able  to  obtain  a  good  estimate  of  the  covariance  matrix  using  fewer,  but  higher  quality,  •<  'ta 
vectors  than  in  the  conventional  approach.  The  DVSO  method  is  applied  in  order  to  And  tf.*  .i4 
best  data  vectors  out  of  a  set  of  N  vectors,  [x,].  This  process  is  repeated  Q  times  to  give 

«r- 

data  vectors  which  are  used  to  calculate  the  covariance  matrix  R: 


The  final  step  applies  the  DVSO  method  to  the  columns  of  the  covariance  matrix  to  And  an 
adapted  weigjit  vector. 

The  number  of  computations  required  by  the  DVSO-COVAR  technique  is  data-dependent; 
the  mean  number  of  calculations  is  estimated  by  Monte-Carlo  simulations.  For  a  single  trial 
evaluation  of  the  projection  matrix,  the  number  of  multiplications  is  [16]: 

0 

1^1 


f 


f 


where 


*0  .SMN 


6.13  Eigenvector-based  projection  technique 

This  technique  is  based  on  an  eigendecomposidon  of  the  covariance  matrix.  The 
covariance  matrix  is  calculated  as  in  Yeh’s  technique  with  equation  (3).  This  Hennitian  matrix 
is  then  decomposed  in  terms  of  its  eigenvectors  as 

R  -  ^  AiViV? 
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where  M  is  the  number  of  data  vectors  used  in  the  formation  of  the  covariance  matrix  and  K  is 
the  number  of  array  elements.  The  second  term,  C„,  is  associated  with  the  calculation  of  the 
1  eigenvalues  and  associated  eigenvectors.  The  total  number  of  multiplications  for  the  evaluation 

of  the  eigenvalues  and  some  associated  eigenvectors  is 

’  ^0v  "  ^er  ■*' 

where  Q.  the  number  of  operations  required  for  the  tridiagonalization,  is  given  by 

''' 


The  number  of  computations  for  the  QL  algorithm  is 

C'ot  -  6K^-3IC, 

The  number  of  operations  required  for  the  power  method,  is 

Cp^  -  ♦  i2Kn^  *  6f,(0  ♦  8n,^ 
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where  n„  is  the  number  of  eigenvectors  to  be  calculated  and  f,(n„’)  is  the  number  of 
computations  required  for  the  full  eigen  solution  of  a  n„  x  n„  matrix. 


6^  Simulation  results  •  projection  techniques 

This  section  compares  SEWE  technique  with  DVSO,  DVSO-COVAR  and  an  eigen 
decomposition  technique  in  terms  of  SNJR  versus  computational  load.  We  performed  simulations 
for  the  three  scenarios  discussed  previously  with  JNR=26.9  dB.  In  DVSO-COVAR,  the 
parameter  N  was  set  to  approximately  (1.5  L),  while  Q  was  varied  from  1  to  20.  Figures  16,  17 
and  18  give  SNJRu,  versus  the  number  of  complex  multiplications  for  DVSO.  DVSO-COVAR. 
SEWE  and  an  eigenvector-based  projection  technique  for  scenarios  1  (L=3),  2(L=5)  and  3  (L=10) 
respectively.  When  L«K  as  in  scenario  1,  the  SEWE  technique  achieves  near  optimum 
performance  (only  .53  dB  below  maximum  value  while  DVSO-COVAR  is  0.2  dB  below)  with 
20  times  fewer  computations  ihan  DVSO-COVAR.  For  larger  L  (scenario  3),  the  SNJRjo  was 
2.72  dB  below  the  optimum  in  SEWE  method  as  opposed  to  .34  dB  for  DVSO-COVAR.  But 
SEWE  required  10  times  fewer  computations  than  DVSO-COVAR.  For  all  the  three  scenarios, 
eigen  decomposition-based  method  achieves  the  optimum  SNJRjo  but  at  the  cost  of  a  huge 
computational  load.  DVSO  gives  sub-optimum  performance  but  with  far  fewer  computations. 

These  results  show  that  SEWE  technique  is  certainly  the  best  choice  for  cases  where 
L«K  (scenarios  1  and  2  for  example).  It  is  also  a  good  choice  for  a  larger  number  of  jammers 
if  one  is  willing  to  trade-off  some  performance  for  a  large  reduction  in  the  number  of 
computations. 
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Complex  Multiplications 

Figure  16  Normalized  SNJR20  versus  Complex  Multiplications  with  Ls3,  K=64  and 


JNR-26.9  dB 


Complex  Multiplications 


Figure  17  Normalized  SNJRjo  versus  Complex  Multiplication  for  L=5,  K=64  and  JNR=26.9 
dB. 
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Complex  Multiplications 


Figure  18  Nonoalized  SNIR^q  versus  Complex  Multiplication  for  La  10,  Ka64  and  JNRs26.9 

dB. 


7.0  Conclusions 

This  report  has  shown  that  the  Yeh  projection  technique  significantly  improves  the 
performance  of  adaptive  beamfonning  provided  a  priori  information  is  available  on  the  number 
of  jammers.  An  improvement  to  the  Yeh  technique  has  been  developed  and  presented;  this 
improved  technique,  called  Subspace  Estimation  witout  Eigenvecton  (SEWE),  does  not  require 
a  priori  knowledge  of  the  jammer  scenario.  The  SEWE  technique  has  been  compared  with  other 
fast  projection  methods  and  also  with  an  eigenvector  decomposition  technique.  The  SEWE 
technique  demonstrated  superior  performance  in  terms  of  Sig^-to-Noise-plus-Jammer  Ratio 
(SNJR)  versus  computational  load,  except  when  a  very  high  SNJR  is  required;  in  this  latter  case, 
DVSO-COVAR  and  eigenvector  projection  perform  somewhat  better.  It  sould  be  noted  however, 
that  the  best  performance  for  SE^l^  for  a  64  element  airay  with  10  jammers  was  a  loss  of  only 
2  dB  with  respect  to  the  optimum.  This  loss  was  considerably  less  for  smaller  number  of 
jammers,  e.g.  3  and  5  jammers  where  the  loss  was  0.2  dB  and  1.2  dB  respectively. 
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